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ABSTRACT 

How much amount and what size of dust are formed in the ejecta of core- 
collapse supernovae (CCSNe) and are injected into the interstellar medium (ISM) 
depend on the type of CCSNe through the thickness of outer envelope. Recently 
Cas A was identified as a Type lib SN (SN lib) that is characterized by the small- 
mass hydrogen envelope. In order to clarify how the amount of dust formed in 
the ejecta and supplied into the ISM depends on the type of CCSNe, we inves- 
tigate the formation of dust grains in the ejecta of a SN lib and their evolution 
in the shocked gas in the SN remnant (SNR) by considering two sets of density 
structures (uniform and power-law profiles) for the circumstellar medium (CSM). 
Based on these calculations, we also simulate the time evolution of thermal emis- 
sion from the shock-heated dust in the SNR and compare the results with the 
observations of Cas A SNR. We find that the total mass of dust formed in the 
ejecta of a SN lib is as large as 0.167 Mq but the average radius of dust is 
smaller than 0.01 /im and is significantly different from those in SNe II-P with 
the massive hydrogen envelope; in the explosion with the small-mass hydrogen 
envelope, the expanding He core undergoes little deceleration, so that the gas 
density in the He core is too low for large-sized grains to form. In addition, the 
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low-mass hydrogen envelope of the SN lib leads to the early arrival of the reverse 
shock at the dust-forming region. If the CSM is more or less spherical, therefore, 
the newly formed small grains would be completely destroyed in the relatively 
dense shocked gas for the CSM hydrogen density of tt-h > 0.1 cm~^ without being 
injected into the ISM. However, the actual CSM is likely to be non-spherical, so 
that a part of dust grains could be ejected into the ISM without being shocked. 
We demonstrate that the temporal evolution of the spectral energy distribution 
(SED) by thermal emission from dust is sensitive to the ambient gas density and 
structure that affects the passage of the reverse shock into the ejecta. Thus, the 
SED evolution well reflects the evolution of dust through erosion by sputtering 
and stochastic heating. For Cas A, we consider the CSM produced by the steady 
mass loss of M ~ 8 x 10~^ Mq yr~^ during the supergiant phase. Then we find 
that the observed infrared SED of Cas A is reasonably reproduced by thermal 
emission from the newly formed dust of 0.08 M©, which consists of the 0.008 Mq 
shock-heated warm dust and 0.072 Mq unshocked cold dust. 

Subject headings: dust, extinction — infrared: ISM — ISM: supernova remnants 
— shock waves — supernovae: general — supernovae: individual (Cassiopeia A) 

1. Introduction 

Core-collapse supernovae (CCSNe) have been considered to be one of important sources 
of interstellar dust. They eject a large amount of heavy elements synthesized during their 
stellar evolution and explosive burning, and form dust grains inside the dense, metal-rich 
He core in the expanding ejecta. Meanwhile, the reverse and forward shocks driven by the 
encounter of the ejecta materials with the surrounding medium propagate into the SN ejecta 
and the interstellar medium (ISM), respectively, producing the hot gas plasma within the 
supernova remnant (SNR). The dust grains formed in the He core, before being dispersed 
into the interstellar space, pass through the shock-heated gas and suffer from erosion and 
destruction by sputtering due to collisions with the energetic ions. Since the efficiency of 
destruction of dust grains heavily depends on their composition and size, what composi- 
tion and what size of dust condense in the ejecta should be clarified for understanding the 
properties and amount of dust finally injected from CCSNe into the ISM. 

Observations of nearby CCSNe have reported increasing evidence for dust formation in 
the ejecta (see Kozasa et al. 2009 for review). Nevertheless, the composition and size of 
newly formed dust have not been specified except for two cases; SN 2004et (Kotak et al. 
2009) and SN 2006jc (Sakon et al. 2009; Smith et al. 2008; Mattila et al. 2008). The masses 
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of dust formed in the SN ejecta are estimated to be less than 10~^ M© from the observations 
by usually assuming that the ejecta is optically thin for thermal emission of dust. However, 
newly formed dust may reside in optically thick clumps (e.g., Miekle et al. 2007). Moreover, 
some species of dust grains can cool down quickly after the formation (Nozawa et al. 2008), 
so that near- to mid-infrared (IR) observations could miss the cold dust. Hence, the derived 
mass of dust could be the lower limit in some cases, and the sophisticated radiative transfer 
calculations (e.g., Sugerman et al. 2006; Ercolano et al. 2007) referring to the observational 
data at UV/optical to far-IR wavelengths would be required to deduce the properties and 
mass of dust formed in the ejecta from observations. 

Theoretical studies have shown that various species of dust grains with the radii of 
0.001-1 /xm and total mass of 0.1-1 Mq can condense in the ejecta of Type 11-P SNe (SNe 
H-P) retaining the massive hydrogen envelope (Kozasa et al. 1989, 1991; Todini & Ferrara 
2001; Nozawa et al. 2003). Adopting the results of dust formation calculations, a semi- 
analytical model (Bianchi & Schneider 2007) and a numerical simulation (Nozawa et al. 
2007) investigated the survival of dust in the SNRs and found that dust mass injected into 
the uniform ISM is ~0.01-0.8 Mq for the ambient hydrogen density of nH,o = 10-0.1 cm~^. 
This mass of dust supplied by SNe II-P is still large enough to explain a huge quantity 
of dust in excess of 10*^ Mq observed for the host galaxies of quasars at redshift z > b 
(Morgan & Edmunds 2003; Maiolino et al. 2006; Dwek et al. 2007), though it is proposed 
that asymptotic giant branch stars can also contribute to dust enrichment even in such an 
early epoch (Valiante et al. 2009). 

Recently Kozasa et al. (2009) have suggested that how much amount and what size of 
dust are formed in the ejecta and are injected into the ISM depend on the type of CCSNe, 
reflecting the thickness of the H-rich envelope. The less massive outer envelope leads to 
extremely high expansion velocity of the He core unless the kinetic energy of explosion is 
very low, and makes the gas density in the expanding ejecta of the He core too low for large- 
sized grains to be produced. Furthermore, the low-mass outer envelope can allow the early 
arrival of the reverse shock at the dust-forming region and result in efficient destruction of 
newly formed dust in the shock-heated gas with high density. Because a significant fraction 
of CCSNe (~0.3) explode as hydrogen/helium-poor CCSNe such as SNe Ib/Ic and SNe lib 
(Prieto et al. 2008; Smartt et al. 2009; Boissier & Prantzos 2009), it is essential to investigate 
how the formation and destruction processes of dust depend on the type of CCSNe, in order 
to reveal the role of CCSNe as sources of interstellar dust. 

It should be mentioned here that dust grains embedded in the postshock gas are eroded 
by sputtering and are heated up by collisions with the hot gas, radiating their thermal 
emission at IR wavelengths. Since the destruction efficiency and heating rate of dust are 
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sensitive to grain size as well as gas density, the IR spectral energy distributions (SEDs) 
by thermal emission from dust in SNRs can be a diagnostic of the size and temperature 
of dust, and therefore provide key information on the nature of dust formed in the ejecta, 
density structure in the ambient medium, and destruction and heating processes of dust. 
In particular, young SNRs that originated from SNe with the low-mass envelope are ideal 
ones to compare the observed IR SEDs with calculated ones; the reverse shock can quickly 
arrive at the He core and heat up newly formed dust before the forward shock sweeps up 
a considerable amount of circumstellar /interstellar dust, so that we can mainly see the IR 
emission from dust grains condensed in the ejecta. In this sense, Cassiopeia A (Cas A) SNR 
is one of the best targets for this study because the type of the Cas A SN is recently identified 
as Type lib by a striking similarity between the optical spectrum of scattered light echo of 
Cas A SN and the time-averaged spectrum of SN lib 1993 J (Krause et al. 2008), and the 
presence of dust freshly formed in the ejecta has been confirmed by IR observations (e.g., 
Rho et al. 2008). 

In this paper we explore the formation of dust in a SN lib and the evolution of the dust 
inside the SNR, with the aim of revealing the composition, size, and mass of dust ejected 
from hydrogen-deficient CCSNe. In addition, the time evolution of thermal emission from 
dust within SNRs is calculated and is compared with the observations of Cas A to discuss 
the composition and amount of dust. In § 2 we perform the calculation of dust formation 
in the ejecta of a SN lib whose progenitor has lost most of the hydrogen envelope during 
the evolution (Nomoto et al. 1993; Podsiadlowski et al. 1993). In § 3 we investigate the 
evolution of the newly formed dust in the gas heated by the reverse and forward shocks and 
its dependence on the ambient gas density and structure. Then, following the calculations of 
dust evolution in SNRs, we present the temporal evolution of IR emission from collisionally 
heated dust in the hot gas in § 4, and compare the results with the observations of Cas A in 
5 5. We summarize our results in 5 6. 



2. Formation of Dust in the Ejecta of SN lib 

In the ejecta of SNe, dust formation takes place within the He core where the densities 
of condensible elements are high enough to advance the formation and growth of seed nu- 
clei through collisions with the relevant gas species. Formation processes of dust grains are 
influenced by the time evolution of the density and temperature of the gas as well as the ele- 
mental composition in the ejecta (Kozasa et al. 1989). In this section, we describe the model 
of the dust formation calculation for a SN lib and present the results of the calculations. 
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2.1. Model of SN lib and Calculation of Dust Formation 

The calculation of dust formation is carried out by applying the theory of non-steady 
state nucleation and grain growth developed by Nozawa et al. (2003). The theory takes 
into account chemical reactions at condensation, considering that the kinetics of nucleation 
and grain growth is controlled by the key species defined as the gas species with the least 
collisional frequency among the reactants (Kozasa & Hasegawa 1987; Hasegawa & Kozasa 
1988). Given the elemental composition of the gas in the ejecta as well as the time evolution 
of the gas density and temperature, the theory allows us to determine the condensation time, 
chemical composition, size distribution, and mass of dust that condenses in the ejecta (see 
Nozawa et al. 2003 for details). 

In the calculation of dust formation, we adopt the model of the SN lib with the zero-age 
main-sequence mass of 18 Mq and metallicity oi Z = 0.02. The ejecta mass is Meje = 2.94 
Mq, the mass of hydrogen envelope is = 0.08 Mq, the kinetic energy of explosion 

is -Ekin = 10^^ ergs, and the synthesized ^^Ni mass is M(^^Ni) = 0.07 Mq. These model 
parameters are invoked to reproduce the light curve of SN 1993J by referring to the SN 
model in Shigeyama et al. (1994). Note that the mass of ejecta in this model is in good 
agreement with the ejecta mass of Gas A (= 2-4 Mq) estimated from the studies of shock 
dynamics (Borkowski et al. 1996; Ghevalier & Oishi 2003; Laming & Hwang 2003) and from 
X-ray observations (Vink et al. 1996; Favata et al. 1997; Willingale et al. 2003). 

In Figure 1 we show the elemental composition in the ejecta of the SN lib model at 
600 days after the explosion, taking the decay of radioactive elements into account. In 
Table 1 we also present the mass of the main refractory species i in the ejecta. We assume 
that the elemental distribution in the ejecta remains an onion-like structure without any 
mixing, which is consistent with the theoretical model of mixing in SN 1993J; Iwamoto et 
al. (1997) have shown with the 2D simulation that the deceleration of the He core by the 
thin H-rich envelope is so weak that the growth of the Rayleigh- Taylor instability is too 
slow to cause substantial mixing in the ejecta, in contrast to SN 1987A which has a thick 
H-rich envelope (e.g., Hachisu et al. 1990, 1992). Although the jet formation in the explosion 
suggested from the observations of Gas A can lead to macroscopic mixing (e.g., Tominaga 
2009), a simple argument of molecular diffusion has shown that mixing at the atomic level 
is extremely ineffective before the condensation of dust grains in the ejecta (Deneault et 
al. 2003). Indeed, the observations of Gas A at various wavelengths have suggested that 
the mixing in the ejecta should not be at atomic level but at the macroscopic level (e.g., 
Willingale et al. 2002; Fesen et al. 2006b; Ennis et al. 2006). From Figure 1, we can see 
the layered structure of elemental composition with the very thin hydrogen envelope at the 
outermost ejecta, which is followed by G-rich (G/0 > 1) He layer, 0-rich layer, Si-S-Fe layer. 
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and Fe-Ni layer inwardly. 

The distributions of gas temperature and density in the ejecta of the SN lib at 300 days 
and 600 days are shown in Figures 2a and 2b, compared with those for the SN II-P model 
of Meje = 17 Mq and Menv = 13.2 Mq with the same Ekin and M(56Ni) as those of the SN 
lib model (Umeda & Nomoto 2002). The gas density scales as with the time t after 
the explosion, and the time evolution of gas temperature is mainly regulated by the energy 
deposited from radioactive elements (i.e. M(^^Ni)). It should be noted that the gas density 
in the He core of the SN lib is about three orders of magnitude lower than that in the SN 
II-P. The reason for this is that, since almost all of the explosion energy are deposited into 
the He core without the thick hydrogen envelope, the resulting expansion velocity of the gas 
in the He core of the SN lib is much higher than that of the SN II-P, as shown in Figure 2 c. 

Based on the hydrodynamic and nucleosynthesis models for the SN lib, we calculate 
dust formation. Dust formation proceeds through the formation and growth of seed nuclei 
as the gas cools down. The number of seed nuclei is evaluated by calculating the non-steady 
state nucleation rate (Eq (5) in Nozawa et al. 2003) being coupled with the time evolution 
of the gas temperature and the number density of the key species. In the calculations, it is 
assumed that without coagulation the seed nuclei grow homogeneously at a rate controlled 
by the collisions of the key species with the sticking probability of unity. We also assume 
that the formation of CO and SiO molecules is complete; all C (Si) atoms are locked in CO 
(SiO) molecules in the case where the number density of C (Si) atoms is lower than that 
of O atoms. Under this assumption, C-bearing grains can condense only in the He layer 
with C/0 > 1. Depletion or rarefaction of the key species due to the grain growth or the 
expansion of the ejecta, respectively, causes the process of nucleation and grain growth to 
cease in less than a few tens of day since the onset of nucleation. 



2.2. Results of Dust Formation Calculation 

The simulation results of dust formation in the SN lib are illustrated in Figures 3a and 
3b, where the mass distribution and average radius of each dust species formed in the ejecta 
are given as a function of enclosed mass Mr, respectively. Without any mixing of elements, 
different species of dust condense according to the elemental composition of the gas in each 
layer: C grains in the He layer, silicates (MgSiOa, Mg2Si04, and Si02) and oxides (AI2O3, 
and MgO) in the 0-rich layer, and FeS and Si grains in Si-S-Fe layer. The typical times (gas 
temperatures) at which dust grains condense in the ejecta are 300-350 days (1400-1700 K) 
for C grains, 350-500 days (1100-1400 K) for silicate and oxide grains, and 600-700 days 
(800-900 K) for Si and FeS grains. In the innermost Fe-Ni layer, neither Fe nor Ni grains 



- 7- 



can be formed because the gas density becomes too low before the gas cools down to the 
temperatures (^800 K) at which these grains can condense. Fe atoms residing in both the 
0-rich and He layers are not abundant enough to be locked in grains as Fe-bearing dust (see 
Fig. 1). Except for Fe grains, the major species of dust formed in the SN lib are the same 
as those for the unmixed cases in SNe II-P with Z = in Nozawa et al. (2003). 

The mass of each dust species formed in the ejecta of SN lib is presented in Table 2, 
together with the ratio of its mass to the total dust mass (see Table 1 for the mass and mass 
fraction of refractory elements locked in dust grains). The total mass of newly formed dust is 
0.167 M0, which is more than two orders of magnitude larger than the dust mass estimated 
from a handful of observations of nearby dust-forming SNe. C grains are the most dominant 
species in mass, with the mass of 0.07 Mq, and the bulk of them is in the innermost He layer 
at Mr = 2.7-2.9 Mq (see Fig. 3a), where the mass fraction of available C atoms locked in C 
grains is /c = 0.4-1. The masses of MgSiOg, Mg2Si04, and SiOa are 0.055 Mq, 0.017 M©, 
and 0.016 Mq, respectively, and these silicate grains account for more than 50 % of the total 
dust mass. The mass distributions of silicate grains are influenced by the abundance of SiO 
molecules (namely Si atoms, see Fig. 1), and fairly concentrate at Mr = 1.7-1.9 Mq with the 
mass fraction of the preexisting SiO molecules locked in silicate grains /si = 0.4-0.8. The 
condensation efficiency, defined as the ratio of the total dust mass to the total metal mass, 
is ~ 0.13 for the SN lib and lies in the range of those estimated for SNe II-P (0.05-0.25, 
Nozawa et al. 2003). 

It should be emphasized here that, as can be seen from Figure 3b, the average radii of 
all grain species formed in the SN lib are below 0.01 /im and are much smaller than those 
(oavc > 0.01 fim) of dust grains formed in SNe II-P (Nozawa et al. 2003). This difference 
is caused by the difference in the gas density within the expanding ejecta of the He core. 
In general, a larger supercooling necessary for dust grains to condense in a gas with lower 
density leads to the formation of a large number of seed nuclei and results in smaller average 
radii even if all available gaseous atoms are used up for the growth of seed nuclei {fj = 1, 
where fj is the mass fraction of key species j locked in a given grain species). As was 
discussed in Kozasa et al. (1989), in the condition that the number density of key species is 
so low that the collision timescale of key species is longer than the expansion timescale of 
the gas, neither a significant number of seed nuclei can condense nor nucleated grains can 
grow to sufficiently large radii via accretion onto the surface of the grains. 
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2.3. Dependence of Grain Radius on SN Parameters 

Here we briefly discuss liow our results on the dust formation are affected by the model 
parameters of SN explosion, in particular focusing on the kinetic energy of the explosion 
-Ekin- A lower E]^i^ leads to a higher He core density, thus making it possible for more dust 
grains with larger radii to condense if the time evolution of the gas temperature could not 
be affected by the explosion energy. At a fixed time ti, the gas density is scaled by the 
expansion velocity Veje as peje(^i) oc v~^^ for a given mass and density structure of the ejecta. 
Since Ve]e oc (-Ekin/Mjjc)^^^ and roughly Oavc oc pl-^ (Nozawa et al. 2003), we find Oavc oc 
1/fcjc oc (MejcZ-Ekin)^^^- Thus, eveu for the extremely low explosion energy of i?kin = 10^° 
ergs, the increase in the average grain radii is as small as a factor of ~ 3. In addition, if we 
consider E'km = (1-4) x 10^^ ergs estimated for Cas A (Willingale et al. 2003; Chevalier & 
Oishi 2003; Laming & Hwang 2003; Young et al. 2006), grain radii only reduce by a factor 
of < 2, compared with those given in Figure 3b. 

In reality, however, the assessed variation of the grain radius is suppressed by the depen- 
dence of the time evolution of gas temperature on -Ekin! with the same amount of ^^Ni in the 
ejecta, the lower (higher) explosion energy can cause the gas temperature to decrease more 
slowly (more quickly) because of more (less) efficient energy deposition in the ejecta with 
the higher (lower) gas density. This leads to the later (earlier) condensation of dust, which 
compensates for the increase (decrease) in the gas density. Consequently, the condensation 
of dust occurs at the similar gas density, and not only the radius but also the mass of dust 
formed in the ejecta is almost the same, regardless of the explosion energy (see Appendix A 
in Nozawa et al. 2003), as long as the range of variation of -Ekin is not extremely large. The 
qualitative analysis indicates that grain radius is insensitive to the explosion energy, and 
that our results on dust formation are robust for the expected variation of model parameters 
for Cas A SN. 

On the other hand, it should be noted that grain radius is greatly affected by the 
structure of the progenitor star resulting from the presence or absence of the massive outer 
envelope; as seen from Figure 2c, the He-core velocity of the SNe lib with the thin hydrogen 
envelope is 5-10 times higher than that of the SN II-P with the thick envelope. The higher 
velocity of the SN lib leads to a few orders of magnitude lower gas density and thus produces 
about one order of magnitude smaller dust grains than the SN II-P. Nozawa et al. (2008) 
have also pointed out the effect of the lack of the outer envelope on the size of newly formed 
dust in Type lb SN 2006jc, applying the SN model by Tominaga et al. (2008). They found 
that too low gas density in the He core, which is caused by the absence of the outer envelope, 
prevents dust grains with the radii larger than 0.01 /im from forming within the He core. 
These results indicate that the radii of dust grains formed in SNe lib and SNe Ib/Ic are to 



- 9 - 



be generally smaller than those in SNe II-P with the thick outer envelope. 

Our results demonstrate that the composition and mass of dust grains formed in the 
ejecta are almost independent of the type of CCSNe and the initial metallicity of the progen- 
itors because the distribution and yield of metals in the He core are similar. On the other 
hand, their radii strongly depend on the type of CCSNe; SNe with the lower-mass envelope 
produce much smaller dust grains (Oavc < 0.01 fim) than SNe II-P with the much more 
massive envelope. As shown in the next section, small sizes of newly formed dust grains will 
determine their subsequent fates within the SNR. 

3. Evolution of Dust in Type lib SNR 

Dust grains formed in the He core initially comove with the gas without suffering from 
any processings. Upon encountering the reverse shock, the gas is quickly decelerated, while 
the dust grains decouple from the gas with high velocities relative to the gas, and intrude 
into the hot plasma region created by the passage of the reverse and forward shocks. These 
shocked dust grains are decelerated by the drag force, eroded by thermal and/or kinetic 
sputtering, and heated up by collisions with the hot gas. The motion, destruction, and heat- 
ing of dust depend on not only the properties of dust but also the density and temperature 
in the postshock gas. In this section, applying the result of the dust formation calculation, 
we examine the evolution of the dust in the hot gas within SNRs and its dependence on the 
gas density and structure in the CSM. 

3.1. Model of Calculations for Destruction and Motion of Dust in SNRs 

The calculations of dust evolution within SNRs follow the method described by Nozawa 
et al. (2007); dust grains are treated as test particles, and the motion and destruction of dust 
are pursued by calculating self-consistently the erosion due to sputtering and deceleration 
due to the gas drag together with the time evolution of gas temperature and density in the 
spherical symmetry shock. We make use of the size distribution and spatial distribution of 
each dust species obtained in the last section as the initial conditions, and adopt the model 
of the SN lib used in § 2 as those for the velocity and density structure in the ejecta. 

For the structure of circumstellar medium (CSM), which affects the evolution of SNRs 
and thus has critical impacts on the evolution of dust, we consider two sets of density profile. 
One is the uniform CSM with the hydrogen number densities of r;,H,o= 0.1, 1, and 10 cm~^, 
assuming that the gas in the CSM has the solar composition and are fully ionized. For the 
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other CSM density profile, we consider the density structure produced by the steady stellar 
wind during the evolution of the progenitor; the density profile is represented in terms of 
hydrogen number density as 

nii(r) = rim ( — ) cm^^ for r < r2 

= 1 cm~^ for r > r2, (1) 

where ri = 1.2 x 10^^ cm is the radius of the freely expanding outermost ejecta at t = 10 
yr after the explosion in the present model, and r2 is the radius such that ni{^i{r2/ri)~'^ = 1 
cm~'^. The coefficient nH,i is given as 

nu,! = niiiri) = - — ^ , (2) 

where M is the mass-loss rate, Vy^ is the stellar wind velocity, /x is the mean molecular weight, 
and rriH and An are the mass and the number abundance of hydrogen atom, respectively. In 
the calculations, we consider three cases of nH,i = 30, 120, and 200 cm~^; for example, the 
case of nH,i= 30 cm~^ with v^, = 10 km s~^ corresponds to M = 2 x 10"^ Mq yr~^, which is 
the mass-loss rate of Cas A's progenitor suggested from the analysis of shock dynamics by 
Chevalier & Oishi (2003). 

The CSM models used in the calculations are listed in Table 3, and the uniform CSM 
and the stellar wind CSM (p oc r~^) are referred to as the model A and B, respectively, 
with the value of ?t,h,o ot ?7,h,i attached for each case. In the calculations, the spatial bin 
for hydrodynamic calculation is linearly divided into 2510 bins, and the size bin of dust is 
logarithmically divided into 80 bins for the range of 1 A to 1 /im. We assume that the gas 
in the CSM has a temperature of 10^ K and that the ejecta collide with the CSM at t = 10 
yr after the SN explosion. We have confirmed that the moderate changes of the CSM gas 
temperature (~10'^-10^ K) do not infiuence the results of calculations. The calculations are 
performed up to t = 10^-10^ yr, depending on the CSM density. 



3.2. Results of Dust Evolution Calculation 

3.2.1. Penetration of the Reverse Shock into the Ejecta 

We first describe how the penetration of the reverse shock depends on the gas density in 
the CSM. Figure 4 exhibits the time evolution of radius and velocity of the reverse shock for 
the uniform CSM models AO.l, Al, and AlO with nH,o = 0.1, 1, and 10 cm~^, respectively. 



- 11 - 



In the upper panel, we also show the trajectory of the boundary between the He core and 
the hydrogen envelope, for a reference, assuming that the boundary is not decelerated even 
after the collision with the reverse shock. 

The interaction of the SN ejecta with the ambient medium generates the forward shock 
running into the CSM to compress the ambient medium. The shocked hot ambient gas 
thrusts the ejecta in turn, driving the reverse shock that propagates into the ejecta but 
initially expands outward in the Eulerian coordinate. For the higher CSM density, the 
forward shock acquires a lower velocity, and the denser shocked gas is produced behind the 
forward shock. Accordingly, the reverse shock travels through the ejecta with a smaller 
outward velocity so that the unshocked ejecta can catch up with the reverse shock at earlier 
times (see Figure 4). Thus, the time at which the reverse shock collides with the He core 
is earlier for the higher gas density in the CSM: tcou = 355 yr, 155 yr, and 75 yr for the 
model AO.l, Al, and AlO, respectively. For the stellar wind CSM, the collision times are 
earlier than for the uniform CSM models because of higher ambient density at the time of 
the ejecta-CSM interaction: tcoii = 63 yr, 48 yr, and 45 yr for the model B30, B120, and 
B200 with riH.i = 30, 120, and 200 cm^^ (see Table 3). 

Note that these collision times for the remnant of Type lib SN (hereafter Type lib SNR) 
are a few tens times earlier than those for Type II-P SNRs, for which tcou = (1.5-10) xlO^ yr 
for the CSM densities of nn.o = 10-0.1 cm~^ (Nozawa et al. 2007). The earlier arrival of the 
reverse shock at the He core in Type lib SNR is due to its much thinner hydrogen envelope 
than that of Type II-P SNR and the resulting higher expansion velocity of the He core (see 
Figure 2 c). As shown below, the thickness of the hydrogen envelope has great effects on the 
evolution of dust in SNRs. 



3.2.2. Destruction of dust for the uniform CSM 

Figure 5 shows the time evolution of gas density {upper panel) and temperature [lower 
panel) within Type lib SNR for the model AlO (uniform CSM with nH,o= 10 cm~^). The 
gas densities behind the reverse and forward shocks, whose positions are depicted by the 
downward- and upward-pointing arrows in the lower panel, respectively, increase by a factor 
of about 4, compared with the unshocked ones. The temperature of the gas swept up by 
both the shocks rises up to >10^ K and remains over 10^ K by t > 10^ yr. Dust grains 
staying in this hot gas are subject to the destruction by sputtering. 

Figure 6 presents the time evolution of positions [upper panel) and radii [lower panel) 
of dust grains with the initial radii of 0.001 /im (for AI2O3 and FeS) and 0.003 /xm (for the 



-VI- 



other grains) for the model AlO, together with the trajectories of the reverse and forward 
shocks in the upper paneL In this model, the reverse shock encounters the He layer at t = 75 
yr, the 0-rich layer at t = 270 yr, and the Si-rich layer at t = 700 yr. Once newly formed 
dust grains intrude into the shocked gas, they are efficiently decelerated by the gas drag and 
are quickly trapped in the hot plasma. These grains are completely destroyed and are not 
injected into the ISM. Actually all newly formed grains are destroyed in the postshock gas 
for most of the models considered in this paper. One of the reasons for this is that the radii 
of newly formed grains are considerably small; since the deceleration rate of dust is inversely 
proportional to the grain radius (e.g., Nozawa et al. 2007), dust grains with the radii of <0.01 
/xm are easily captured by the shock-heated gas and are destroyed by thermal sputtering. 
Another reason is that the reverse shock can reach the He core more than 20 times earlier for 
SN lib than for SNe II-P as mentioned above; although the initial gas density in the ejecta 
of the SN lib is three orders of magnitude lower than those in SNe II-P (see Fig. 2b), at 
such an early time, the gas density in the He core of the SN lib is still higher by a factor of 
~10 than those in SNe II-P at several thousands years. This leads to the frequent collisions 
between dust and gas and causes the efficient deceleration and destruction of dust. 

In Figure 7a we compare the time evolution of the total dust mass for the uniform GSM 
with nH,o= 0.1, 1, and 10 cm~^. We can see that dust grains are destroyed more efficiently 
and quickly for higher GSM density. This is because the increase of the gas density in 
the GSM leads to the earlier arrival of the reverse shock at the dust-forming region, which 
keeps the gas density in the He core high enough for dust grains to be eroded by sputtering 
effectively. For nH,o= 0.1 cm~'^, G grains of ~10~^ Mq can survive at t = 10^ yr, while all 
dust grains are destroyed within t = 4 x 10^ yr for nH,o= 1 cm"'^ and t = 5 x 10'^ yr for nH,o= 
10 cm"^. 



3.2.3. Destruction of dust for the stellar wind CSM 

Figure 76 shows the evolution of the total dust mass in SNRs for the stellar wind GSM 
(p oc r~^) models. For these models, dust grains are completely destroyed before t = 10^ yr. 
However, the evolution of dust in SNRs shows the different behavior as shown in Figures 8a 
and 8b, where the time evolutions of mass of each dust species for the model B30 and B120 
are presented. For the model B30 with nH,i= 30 cm~^, a small amount of G grains survive 
after t = 10^ yr, whereas the other grain species are completely destroyed before t ~ 10^ 
yr. In the case of the model B120 with the GSM gas density four times higher than that for 
B30, G grains are destroyed quickly after the collision with the reverse shock, whereas the 
masses of some dust species such as Si02 and Si grains decrease more slowly than those for 
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B30. 

In order to explain these different behaviors of dust evolution, we show the time evolu- 
tions of density structure within SNRs for B30 and B120 in Figure 9a and 9b, respectively. 
Since the deceleration and erosion rates are dependent on the number density of ion, the 
time evolutions of the total ion number density, n, are given in the lower panel in Figure 9, 
along with the positions of C, Si02, and Si grains with the initial radii of 0.003 /im. We also 
present the time evolutions of positions and radii of these C, Si02, and Si grains in Figure 
10. 

For the model B30 (Figs. 9a and 10a), the SNR expands into the CSM with the density 
profile of p oc during only 340 yr, and then propagates into the ISM with the constant 
density of nn = 1 cm~^. Thus, one might guess that the evolution of dust in the SNR may 
be quite similar to that for the model Al with the uniform CSM of nH,o = 1 cm~^. However, 
as presented in Figure 7, dust grains in the ejecta are completely destroyed about 10 times 
earlier for B30 than for Al, because the gas density behind the reverse shock for B30 is about 
8 times higher than for Al at t < 340 yr due to the fast passage of the reverse shock into 
the ejecta (see Table 3). In addition, the density of gas in the region swept up by the reverse 
shock at the initial phase remains to be quite high {n > 0.3 cm""^) even after t = 3000 yr (see 
Fig. 9a; hereafter, such a region is referred to as the dense gas region). As can be seen in Fig. 
10a, dust grains are trapped in the dense gas region shortly after their collisions with the 
reverse shock and are completely destroyed by thermal sputtering. C grains formed in the 
inner C-rich region of the He layer collide with the reverse shock at an early epoch [t ~ 330 
yr) and are eroded by sputtering in the denser gas than those for Si02 and Si grains formed 
inside the 0-rich layer. However, C grains can survive for a relatively long period, compared 
with Si02 and Si grains (see Fig. 8a), because the erosion rate of C grains is smaller than 
the other grain species by a factor of about 5 (Nozawa et al. 2006). 

In the case of the model B120 (Figs. 9b and 10b), the reverse shock more quickly pen- 
etrates into the ejecta than B30, which results in the extended shocked region between the 
reverse shock front and the dense gas region (see Fig. 9b). C grains collide with the reverse 
shock at t 2^ 140 yr and are immediately destroyed in the shocked gas with (330/140)'^ ~ 13 
times higher density than that for B30. Si02 and Si grains, which encounter the reverse 
shock at t ~ 400 yr and t ~ 700 yr, respectively, intrude into the 30-40 times denser post- 
shock gas than the cases for B30. However, the gas density is not high enough for the Si02 
and Si grains to be trapped and completely destroyed. Instead, upon collisions with the re- 
verse shock, these grains are greatly decelerated by the gas behind the reverse shock, which 
prevents them from attaining the dense gas region in a short time; Si02 grains are finally 
trapped in the extended shocked region and are destroyed without arriving at the dense gas 
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region. Si grains are eroded slowly passing through the extended shocked region until they 
can reach the dense gas region at 10^ yr. 

For the model B200, the reverse shock can propagate into the ejecta even earlier than 
B120. Hence, the gas density behind the reverse shock is so high that all grain species are 
trapped and are quickly destroyed once they are swept up by the reverse shock. These results 
demonstrate that the evolution of dust within SNRs strongly depends on the gas density and 
structure in the CSM, as well as on the composition of dust grains and their initial position 
in the ejecta. 

In Table 3 we summarize the total mass Msurv of surviving dust in Type lib SNRs at 
t = 10^ yr for the CSM models considered in this paper. For all CSM models except for the 
model AO.l, dust grains formed in the ejecta are completely destroyed before t = 10^ yr. 
Even for the model AO.l, the destruction efficiency t], defined as the ratio of Mgu^v to the 
initial dust mass Md, is less than 10~^, in contrast to the case of SNe II-P with the thick 
hydrogen envelope, for which 2-80 % of dust grains can survive for the ISM with nH,o = 

10- 0.1 cm~^ (Nozawa et al. 2007). As shown above, the difference in the evolution of dust 
between SNe lib and SNe II-P SNRs is caused by the thickness of outer envelope; for SNe 

11- P retaining the massive hydrogen envelope, the radii of newly formed grains are larger 
than 0.01 fim, and the arrival time of the reverse shock at the dust-forming region inside the 
He core is retarded, both of which act against the destruction of dust grains by sputtering. 
Therefore, how much amount of dust is injected into the ISM strongly depends on the type 
of SNe through the thickness of outer envelope, even if a large amount of dust grains can 
form in the ejecta of CCSNe. From the results of the present calculations, we conclude 
that low-mass envelope SNe such as SNe lib and SNe Ib/c cannot be main contributors of 
interstellar dust as long as the gas density in the CSM is not too low (rin ^0.1 cm~^). 



4. IR Thermal Emission from dust in Type lib SNRs 

Dust grains bathed in the shocked gas inside SNRs are heated up by collisions with 
the hot gas and emit their thermal radiation at IR wavelengths. If the coUisional heating 
rate and radiative cooling rate of dust are balanced, the dust can acquire the equilibrium 
temperature. However, in a rarefied hot plasma such as SNRs, small-sized dust grains 
undergo stochastic heating, which infiuences the emissivity and the resulting SED (e.g., 
Dwek 1986). In this section, based on the results of simulation of dust evolution presented 
in § 3 and taking account into stochastic heating, we calculate the SED by thermal emission 
from shock- heated dust within SNRs. The SNRs are placed at a distance of 3.4 kpc in order 
to compare with the observations of Cas A SNR in the next section. 
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4.1. Calculation of stochastic heating 

The heating and thermal emission of dust are calculated by using the output of the dust 
evolution simulations: density and temperature of the gas as well as size, velocity relative 
to the gas, and number density of each dust species in each spatial bin as a function of 
time. We employ a Monte-Carlo method to calculate the temperature distribution functions 
of dust grains by stochastic heating. In the calculations, we consider only collisions with 
electrons for the heating of dust; under the present situation in which dust grains are quickly 
trapped in the hot gas, collisions with ions are much less frequent, and the heating of dust 
is dominated by collisions with electrons. 

The fraction of incident energy that is deposited into a grain is calculated following the 
prescription by Dwek (1987a), and the stopping range of an electron for each grain species is 
calculated by the approximation formula which we derived based on those by Tabata et al. 
(1972) and Iskef et al. (1983). The heat capacity of each grain species is taken from Takeuchi 
et al. (2005), and the optical constants are taken from the following sources: C (Edo 1983), 
AI2O3 (Toon et al. 1976 for A < 8 /im and Begemann et al. 1997 for A > 8 /im), MgSiOs 
(Dorschner et al. 1995, but the optical constants of Mg2Si04 are used for A < 0.3 fim, see 
Hirashita et al. 2008), Mg2Si04 (Semenov et al. 2003), SiOg (Philipp 1985), MgO (Roessler 
& Huffman 1991), FeS (Semenov et al. 2003), Si (Filler 1985). We re-bin the size of dust 
into 12 bins from 1 A to 0.1 /im and divide the shocked region into 20 spatial bins weighted 
by the mass of the gas, because the calculation is too time-consuming. We have verified that 
even if the number of re-binned size (spatial) bin increases to 20 (50), the resulting SEDs 
are virtually the same. 



4.2. Time evolution of thermal emission from dust 

Figures 11-14 show the time evolution of SEDs by thermal emission from shocked dust 
for different CSM models every 200 yr in the range from 200 yr to 1800 yr after the explosion. 
In these figures, the solid curves display the SEDs taking into account stochastic heating, 
and the dashed curves the SEDs calculated by assuming the equilibrium without stochastic 
heating. We can see that the SEDs taking account of the stochastic heating are strongly 
different from those with the equilibrium temperatures, especially at wavelengths shorter 
than A = 20 /im. This implies that the effect of stochastic heating is very important in 
modelling the IR SED by dust emission within SNRs, as demonstrated by Dwek (1986). 

For the model Al with the uniform gas density (Fig. 11), the reverse shock hits the He 
core at t = 155 yr after the explosion. The mass of C grains formed in the outer C-rich 
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region of the He layer is too small (<5 x 10~^ ^0) to emit the thermal emission significantly 
during t = 200-400 yr. At t > 400 yr, thermal emission from a large amount of C grains 
formed in the inner C-rich part of the He layer always dominates the SED, because the gas 
density in the SNR is not high enough to destroy C grains efficiently. The reverse shock 
reaches the 0-rich layer at t ^ 600 yr, but the emission features of silicate and oxide grains 
are inconspicuous, being overwhelmed by thermal emission from C grains. 

For the higher CSM density of the model AlO (Fig. 12), C grains are more efficiently 
destroyed by sputtering than those in the model Al, and the emission features from other 
dust species get prominent with time. The reverse shock encounters the 0-rich layer at 
t = 270 yr, and the emission feature of Mg-silicates (MgSiOs and Mg2Si04) around A = 10 
/xm begins to be noticeable at t = 600 yr. The feature of Si02 at A = 20 /im is outstanding 
after t = 800 yr, and the emission features from Si and FeS grains appear after t = 1200 yr. 
Because of the efficient destruction of dust, the thermal emission continues to decrease with 
time. 

In the case of the model B30 with the stellar wind's density profile (Fig. 13), the reverse 
shock collides with the He core at t = 63 yr, which is earlier than in the model AlO. Thermal 
emission from C grains dominates the SED by t ~ 1000 yr, and after then the contributions 
from Mg-silicates become recognizable. However, stochastic heating of small C grains still 
contributes to the SED at A < 10 fim, because some C grains can survive even after t = 1000 
yr owing to the low gas density (tt-h = 1 cm~^) at a large distance in comparison with the 
model AlO. We can see that the flux densities signiflcantly increase between 1200 yr and 
1400 yr in this model. This is due to the fact that during this time interval, the reverse 
shock has swept up the ejecta region where a large amount of silicate grains are formed 
{Mr ^1.8 Mq, see Fig. 3a). The dust grains swept up by the shocks continue to be eroded 
by sputtering, while the reverse shock continues to intrude into the inner He core and to 
supply newly formed dust grains to the shocked gas. Thus, the time evolution of SEDs is 
determined by the balance between the amount of dust destroyed in the hot gas and the 
amount of dust swept up by the reverse shock. 

For the model B120 (Fig. 14), the reverse shock can quickly propagate into the ejecta, 
and arrives at the He core, 0-rich layer. Si-rich layers at t = 48 yr, 150 yr, and 560 yr, 
respectively. Thus, we can see the quick evolution of the SED characterized by the features 
of various grains species. Since the gas density is signiflcantly high, especially at an early 
phase, C grains are quickly destroyed as can be seen from the evolution of SED at A < 8 
fj,m, and the excess by stochastic heating of small C grains disappears at t > 600 yr. 

The results of calculations clearly demonstrate that the time evolutions of IR emission 
from dust heated up in SNRs are sensitive to the density of gas as well as the density proflle 
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in the surrounding medium. The ambient gas density controls the passage of the reverse 
shock into the ejecta and affects the destruction by sputtering and stochastic heating of dust 
in the hot gas. Hence, the resulting IR SED well reflects the evolution of dust in SNRs and 
also holds key information on the density structure in the ambient medium. In other words, 
the temporal evolution of SED from dust in shock-heated gas can be a very useful probe for 
investigating the properties of dust in SNRs as well as the density structure in the CSM and 
thus the mass-loss history of the progenitor star during the late phase of its evolution. 

5. Applications to the Cas A SNR 

In this section, we apply the results of calculation for dust evolution and the resulting 
SEDs presented in § 3 and § 4 to Cas A SNR, to give an insight into the composition and 
mass of dust residing in Cas A SNR and to get an implication on the mass-loss rate of the 
progenitor. After briefly summarizing the observed characteristics of Cas A relevant to this 
study, we present the comparison of the calculated SEDs with the observed one. 

5.1. Cas A SNR 

Cas A is one of the best studied Galactic SNRs at a broad range of wavelengths spanning 
from radio to 7-ray. Its estimated age is ~330 years (Fesen et al. 2006a; Thorstensen et al. 
2001), and the distance to Cas A is determined to be 3.4 kpc from the proper motion of 
fast moving knots (FMKs, Reed et al. 1995). The X-ray observations with the Chandra have 
resolved the radii of the reverse and forward shocks to be 95"±10" and 153"±12", respectively 
(Gotthelf et al. 2001), which correspond to R^^^ = 1.57±0.17 pc and = 2.52 ±0.19 pc 
at a distance of 3.4 kpc. 

Recently Krause et al. (2008) obtained an optical spectrum of Cas A SN around the 
maximum brightness via the scattered light echo, and concluded that the Cas A SN was of 
Type lib, originated from a supergiant that had lost most of its hydrogen envelope prior to 
the explosion as has been suggested by the close binary scenario for SNe lib (e.g., Nomoto 
et al. 1995). The numerical simulations (Chevalier & Oishi 2003; Young et al. 2006) and the 
discovery of FMKs containing some hydrogen (Fesen & Becker 1991; Fesen 2001) also seem 
to strongly support this conclusion. 

The presence of dust formed in the ejecta has been confirmed by the IR observations 
with the IRAS (Dwek et al. 1987b), the ISO (Lagage et al. 1996), and the Spitzer (Hines et 
al. 2004). The IR spectra of FMKs taken by the ISO allowed us to examine the composition 
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of dust, leading to the conclusion that the main species of newly formed dust in Cas A are 
silicate and/or oxide grains (Arendt et al. 1999; Douvion et al. 2001). 

There has been a continuing controversy about the amount of dust residing in Cas A 
SNR; the dust mass derived from the observations described above was smaller than 3 x 10"'^ 
Mq. Rho et al. (2008) estimated the mass of dust to be 0.02-0.054 Mq from the spectral 
mapping covering the entire region of Cas A by the Spitzer. Based on the submillimeter 
(submm) observation with the SCUBA, Dunne et al. (2003) have reported the detection of 
2-4 Mq of cold dust in Cas A. Dwek (2004) argued that this cold dust emission is ascribed 
to metallic needles of <10~^ Mq formed in the ejecta. On the other hand, the MIPS 160 /xm 
map of Cas A, complemented by CO and OH line observations, have demonstrated that the 
detected submm emission is originated from interstellar dust in the foreground molecular 
clouds (Krause et al. 2004), and the mass of cold dust associated with Cas A is estimated 
at most to be 0.2 Mq (Krause et al. 2004) and 1.5 Mq (Wilson & Bartha 2005). 

Although the existence of a peculiar species of cold dust with ~1 Mq is also proposed 
from submm polarization (Dunne et al. 2009), we do not here consider the measurements 
at wavelengths longer than 160 /xm. In what follows, we compare the calculated SEDs with 
the synchrotron-subtracted IR flux densities at A = 8-100 fim taken from the photometric 
observations tabulated in Hines et al. (2004). 

5.2. Comparison with the observations of Cas A 

5.2.1. IR emission from shock-heated newly formed dust 

We first compare the results for the uniform CSM with the observations of Cas A. For 
"h,o = 1 cm~'^, IR emission is dominated by stochastically heated small C grains throughout 
the age of the remnant, as shown in § 4.2. This causes the excess emission over the observed 
SED at A < 15 /xm until t = 1000 yr, and the calculated SED cannot reproduce the obser- 
vational data. The fact that the C grain-dominated IR SED does not match with the IR 
measurements of Cas A is also pointed out by Bianchi & Schneider (2007). 

For the case of r2H,o = 10 cm~^, the calculated SED at t = 800 yr is in good agreement 
with the IR observations of Cas A (see Fig. 12). The mass of shock-heated dust at t = 800 
yr is 0.048 Mq, which is consistent with the dust mass of 0.02-0.054 Mq in Cas A estimated 
from the observations with the Spitzer (Rho et al. 2008). The comparison strongly suggests 
that the dominant species of dust responsible for the observed SED are silicate grains and 
that the reverse shock has already reached the inner 0-rich layer in the ejecta. However, 
the remnant age of 800 yr is much older than the age of Cas A, and the SED at t = 330 yr 
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does not agree with the observations. We also found that any uniform CSM models with gas 
density higher than rin.o = 10 cm~^ cannot reproduce the observed SED, partly because the 
dust grains are quickly destroyed and partly because the temperature of dust is too high. 
Thus, we conclude that thermal emission of dust grains heated up in SNRs expanding into 
the uniform CSM cannot explain the observed IR SED of Cas A. 

Next we move to the comparison of the results for the CSM density structure of p cx r~^. 
Figures 15a-15c present the calculated SEDs at t = 330 yr for n-n i = 30, 120, and 200 cm~^. 
In these figures, we also show the contribution of each dust species to the total flux densities. 
Note that the SEDs calculated without stochastic heating are not shown in these figures, 
since they can never reproduce the observations of Cas A for any CSM models. 

The case of nH,i = 30 cm~^ corresponds to M = 2 x 10~^ Mq yr~^ for v^, = 10 km 
s~^, and the density structure has been suggested by Chevalier & Oishi (2003) to reproduce 
the observed radii of the reverse and forward shocks derived from the X-ray observations 
(Gotthelf et al. 2001). Our simulation also approximately reproduces the observed shock 
radii for Cas A: 1.5 pc for the reverse shock and 2.2 pc for the forward shock. However, in 
this case, the reverse shock has not yet reached the 0-rich layer at t = 330 yr (see Table 3). 
The calculated flux densities from the shock-heated C grains of ~0.05 Mq are much higher 
than the observations at the whole mid-IR wavelengths. 

With riH,! = 120 cm~^, the calculated SED reasonably reproduces the observed SED 
£1^^ ^ 30 fim, although there is some discrepancy at A 12 fim. The SED also predicts 
the considerably lower flux densities at A > 50 fim, which can be recovered by considering 
the unshocked dust as will be discussed later. The dust species most contributing to the 
SED is MgSiOs, which is one of the major grain species in Cas A presumed from the IR 
spectra (Douvion et al. 2001; Rho et al. 2008). Mg2Si04, Si02, and MgO grains more or 
less contribute to the SED, and C grains partly contribute at A < 15 /im as well. The total 
mass of shocked dust grains is 0.008 Mq, which is smaller than the estimate by Rho et al. 
(2008) but is larger than those by the other IR observations. 

For TT-H,! = 200 cm~^, the calculated SED does not match with the observations, es- 
pecially at A = 20-60 fim. In this case, the reverse shock invades deeper layers in the 
ejecta and sweeps up larger mass (0.028 Mq) of silicate grains, compared with the case of 
nH,i = 120 cm~^. Therefore, the derived SED produces the much more flux densities than 
the observations. We also calculate the thermal emission from shocked dust for the denser 
CSM than rin,! = 200 cm~^, but the outcomes suffer from much less flux densities at A > 30 
/im because of the efficient destruction of dust and the relatively high temperature of dust. 
In addition, in such dense cases, the shock radii are much smaller than the observations. 
Hence, we conceive that the cases of nH,i ^ 200 cm~'^ are unlikely as the density structure 
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around Cas A. 

In summary, our results indicate that the stellar wind's density profile with ?t,h,i — 120 
cm~^ is most favorable to account for the observed IR emission at the age of Cas A. The 
corresponding mass- loss rate 8 x 10^^ Mq yr ^ is at most a factor 4 higher than the mass-loss 
rates estimated for SNe lib: (3-4)xl0~^ Mq yr~^ (Suzuki & Nomoto 1995; Fransson et al. 
1996) deduced for SN 1993J (at r ~ 1 x 10^^ cm) and 2 x 10"^ Mq yr"^ for Cas A (Chevalier 
& Oishi 2003). In view of the uncertainties involved in the density gradient and dumpiness 
in the CSM, these estimated mass loss rates can be regarded as being consistent. On the 
other hand, the result for tt-h,! = 120 cm~^ predicts quite smaller shock radii than those for 
Cas A (see Table 3). The discrepancy in the shock radii as well as the mass-loss rate could 
be considered to be partly caused by the difference in the explosion energy and the thickness 
of the hydrogen envelope of the SN model used in this study. 

It should be emphasized that the IR emission from newly formed dust grains that are 
heated up in the hot plasma inside SNRs can reproduce the mid-IR (A = 8-30 /im) SED of 
Cas A without any fine-tuning. On the other hand, the obtained results for rin i = 120 cm"'^ 
significantly underestimate the flux densities at longer wavelengths (A > 50 fim). Although 
the thermal emission from cooler dust in the dense optically emitting knots is expected 
to mitigate this difference, modelling of such a dust emission may need a sophisticated 
radiative transfer simulation and is beyond the scope of this paper. Hereafter, we discuss 
the contributions from the following two potential sources to the observed SED at A > 50 
/xm: (1) the circumstellar (CS) dust grains that was formed in the mass-loss wind during 
the evolution of the progenitor star and then were swept up by the forward shock, and (2) 
the newly formed dust that were formed in the ejecta but have not been yet swept up by 
the reverse shock. 



5.2.2. Contribution from preexisting CS dust 

In order to examine the contribution from the preexisting CS dust, we simulate the 
destruction and heating of the shocked CS dust for the density structure of the CSM cor- 
responding to the mass-loss rate of 8 x 10"^ Mq yr~^. The composition of the CS dust is 
considered to be either carbonaceous or silicate depending on the C/0 ratio in the stellar 
wind. Carbonaceous grains condense in C/0 > 1, which is realized in the winds of carbon- 
rich Wolf-Rayet stars. Silicate grains are expected to form in the wind of a supergiant, where 
the abundance ratio of heavy elements is nearly solar with C/0 < 1. Although a supergiant 
is a more likely progenitor of Cas A (§ 5.1), we consider here, as the dust composition in 
the CSM, both amorphous carbon and astronomical silicate, for which the optical constants 
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are taken from Edo (1983) and Draine & Lee (1984), respectively. For the size distribution 
of dust, we adopt the MRN size distribution {n{a)da oc a~^'^da, Mathis et ah 1977) with 
the upper and lower limits of a = 0.001 fim and 0.5 fim, respectively. The dust-to-gas mass 
ratio in the CSM is assumed to be constant and is taken as a free parameter. 

Figure 16 illustrates the contribution from CS carbon (Fig. 16a) and astronomical sili- 
cates (Fig. 16b) for the case of nH,i = 120 cm~^, where the ejecta-dust and CS dust compo- 
nents are depicted by the dotted and dashed lines, respectively. The dust-to-gas mass ratios 
adopted here are 1.6 x 10"'' for carbon and 2.8 x 10~^ for silicate because higher dust-to-gas 
mass ratios cause the excess emission over the observed SED at A < 50 fim. At t = 330 yr, 
the masses of shock-heated CS carbon and silicate grains are 1.4 x 10^^ Mq and 1.8 x 10~^ 
Mq, respectively, and both grain species have typical temperature of ~100 K. Thus, we can 
see that both grain species can only contribute to the SED at A < 50 fim with the peaks 
around A = 30 fim, and cannot compensate the flux densities at A > 50 fim. This means 
that the observed SED cannot be improved by the contributions from CS grains even if 
we consider any mixture of these CS grains and dust-to-gas mass ratios different from those 
considered here. Therefore, we suppose that the IR emissions from shocked CS grains should 
not significantly contribute to the IR SED of Cas A and that a small content of CS dust may 
be due to low condensation efficiency of dust in the mass- loss wind and/or due to partial 
destruction of CS dust by the strong radiation associated with the shock breakout at the 
explosion (e.g., Blinnikov et al. 1998). We emphasize that the thermal emission from CS 
dust alone can never reproduce the observations of Cas A for any models considered in this 
paper. 



5.2.3. Contribution from unshocked ejecta-dust 

As shown in § 5.2.1, in the case of riH.i = 120 cm~'^, only 0.008 Mq of newly formed dust 
grains is heated in the shocked gas at t = 330 yr, and 0.072 Mq of dust grains remains cold 
without being swept up by the reverse shock (see Table 3). These unshocked dust grains 
are expected to be heated up by absorbing the X-ray and UV/optical photons generated 
from the shock-heated plasma and dense knots, and they may radiate their thermal emission 
with the higher temperature than that (~20 K) in the diffuse ISM and contribute to the 
SED at IR wavelengths. In fact, the presence of cold dust in the unshocked ejecta of Cas 
A SNR has been suggested from the association of a "featureless dust" component with the 
Si II emission from the central part of the remnant (Rho et al. 2008). Since the dedicated 
calculations of radiative transport at X-ray to optical wavelengths is beyond the scope of 
this paper, we evaluate the contribution of the unshocked dust to the IR SED, assuming the 
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same temperature for all unshocked grains. 

Figure 17 shows the contribution from unshocked newly formed dust, together with 
the contribution of shock-heated warm dust at 330 yr (Fig. 17a) and 340 yr (Fig. 17b) for 
riH,! = 120 cm~^. The age of 340 yr is the expected age of Cas A in the case where the 
forward shock has not been decelerated (Thorstensen et al. 2001; Fesen et al. 2006a). As 
can be seen from the figures, if the unshocked dust grains are heated up to a temperature 
around 40 K, their thermal emission can significantly contribute to flux densities at A > 50 
/xm and can successfully reproduce the observations of Cas A. Thus, the IR SED observed 
for Cas A can be naturally explained by the shocked warm component and unshocked cold 
component of newly formed dust obtained by a series of our calculations. We also mention 
that the calculated IR emission from shocked warm dust at t = 340 yr slightly increases 
from that at t = 330 yr, giving a better fit at A ~ 20 /im. This indicates that the IR flux 
from Cas A can change during a short time interval of ~10 yr. 

Finally, it should be pointed out here the foUowings: We show that the IR observations 
of Cas A are successfully reproduced by 0.008 Mq shocked dust and 0.072 Mq unshocked 
dust, considering the destruction and heating of newly formed dust in the SNR for n-n i = 120 
cm~^. Although the derived dust mass in Cas A SNR is superficially large enough to explain 
a vast amount of dust observed in the host galaxies of high-redshift quasars at z > 5, these 
dust grains would be completely destroyed in 2 x 10^ yr and could not be supplied into the 
ISM, as demonstrated in § 3.2, if the CSM would be more or less spherical. We also should 
note that if the CSM would not be spherical, a part of dust grains in the unshocked ejecta 
could be ejected into the ISM. 

6. Conclusion 

We investigate the formation of dust in the ejecta of a hydrogen-poor SN lib, adopting 
the SN model with E'kin = 10^^ ergs, Meje = 2.94 Mq, and Menv = 0.08 Mq. Then we explore 
the evolution of newly formed dust in the hot gas inside the SNR by considering the two sets 
of density profile for the spherical CSM. Based on the results of dust evolution in SNRs, we 
also calculate the time evolution of thermal emission from the shocked dust within the SNR 
and compare the results with the observations of Cas A. The results of our calculations are 
summarized as follows. 

(1) In the unmixed ejecta of the SN lib, various dust species can condense depending 
on the elemental composition in each layer at 300-700 days after the explosion. The total 
mass of newly formed dust is as large as 0.167 Mq, and the composition and mass of dust 
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formed in the SN lib are similar to those in SNe II-P. However, the mass of the hydrogen 
envelope of the SN lib is too small to significantly decelerate the expansion of the He core. 
Thus, the gas density in the He core is too low for large-sized grains to be produced, and 
the resulting average radius of dust is confined to be smaller than 0.01 fim. We propose that 
the size of dust formed in the SNe with the low-mass envelope such as SNe lib as well as the 
envelope-stripped SNe Ib/Ic is much smaller than that formed in SNe II-P with the massive 
hydrogen envelope. 

(2) The evolution of dust within SNRs is strongly affected by the mass of the outer 
envelope and the gas density in the CSM. The less massive outer envelope and/or the denser 
CSM lead to the earlier arrival of the reverse shock at the He core, and causes the quicker 
and more efficient destruction of newly formed dust in the denser hot gas. In the Type lib 
SNR, as a result of the small size of dust and the early propagation of the reverse shock 
into the dust-forming region, all newly formed dust grains are destroyed in the shocked gas 
without being injected into the ISM for the ambient gas density of rin > 0.1 cm~^. If the 
CSM is more or less spherical, our results imply that low-mass envelope SNe do not play a 
major role in the dust enrichment in the ISM. 

(3) The time evolution of IR SEDs by thermal emission from shocked dust in SNRs 
well reflects the evolution of dust through the erosion by sputtering and stochastic heating. 
The destruction and heating processes of dust in the hot gas are heavily affected by the 
ambient gas density that affects the passage of the reverse shock into the ejecta. Thus, the 
comparison of theoretically modelled SEDs and observed IR SEDs can be used to diagnose 
the size and mass of dust in SNRs and the density structure surrounding the SNe. These 
in turn provide us with some hints on the formation process of dust in the ejecta and the 
mass-loss activity of the progenitor stars during the late evolution. 

(4) The IR SED observed for Cas A can be reasonably reproduced by the shocked warm 
dust with Md = 0.008 Mq and unshocked cold dust with Ma = 0.072 Mq, both of which 
were newly formed in the ejecta. The CSM density required to reproduce the observed SED 
corresponds to the steady stellar wind of a supergiant with the mass-loss rate of 8 x 10~^ 
Mq yr^^. This mass-loss rate is consistent with those estimated for the progenitors of SN 
1993J and Cas A within a factor of 2-4. It should be emphasized that these conclusions 
are the natural outcomes derived by the self-consistent treatments of formation, destruction, 
and heating of dust given in this study. 

The result of the dust formation calculation shows that the total mass of dust formed 
in the ejecta is 0.167 Mq, which is much larger than the mass estimate of dust deduced from 
the observations of nearby dust-forming SNe. From the comparison of calculated SEDs with 
the observed SED, we predict that in Cas A SNR about a half of newly formed dust grains 
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has already been swept up by the reverse shock and most of them had been destroyed, while 
the other half has not yet been swept up by the reverse shock. The mass of this unshocked 
dust (0.072 Mq) is in good agreement with the upper mass limit of cold dust (0.2 Mq) in Cas 
A estimated by Krause et al. (2004). On the other hand, our calculation demonstrates that 
all dust grains formed in Cas A might be completely destroyed later on and might not be 
injected into the ISM, if the CSM would be more or less spherical (see below). Thus, the SN 
lib that formed Cas A might not be a representative of SNe as major sources of interstellar 
dust, but Cas A could have remained as a unique and ideal object for investigating formation 
and evolution processes of dust in SNe. 

We should bear in mind that all calculations here are performed for the spherically 
symmetric forward and reverse shocks propagating into the homogeneous CSM and stratified 
ejecta, respectively. In the binary scenario for SNe lib (Nomoto et al. 1995, 1996), the 
formation of aspherical CSM is predicted to be as follows. The close binary consists of 
a massive SN lib progenitor and a small-mass companion star. In the common envelope 
phase, the companion star spirals into the envelope of the SN lib progenitor, which causes 
extensive mass loss from the H-rich envelope of the progenitor in the direction of equatorial 
plane. The resultant CSM might have a ring-like structure as seen in SN 1987A, a torus, 
or a thick disk. In reality, the CSM around SN 1993J is considered to be more like torus 
shape (e.g., Fransson et al. 2005). Then since the ejecta of polar direction cannot strongly 
interact with the CSM (e.g., Maeda et al. 2002), the ejecta of SN 1993 J can be composed of 
the reverse- shocked ejecta and the unshocked ejecta. Thus, if CSM is spherical, dust can be 
totally destroyed. If CSM is non-spherical like SN 1993J, unshocked dust can partially be 
ejected. 

In addition, the IR and X-ray observations of Cas A have shown that the ejecta of 
Cas A are very clumpy and asymmetric (Willingale et al. 2002; Ennis et al. 2006; Smith et 
al. 2009) with the inverted compositions as represented by Fe-rich materials overtaking the 
outer layer (e.g., Hughes et al. 2000). This indicates that some dust species may encounter 
the reverse shock at different times and different gas densities than suggested in this paper. 
Furthermore, if newly formed grains would reside in dense clumps, it might be possible that 
a part of them can survive the destruction and be injected into the ISM. Although the mass 
of shocked and visible material is believed to roughly account for all the mass ejected in the 
explosion, the IR spectral mapping observations of the entire region of Cas A SNR suggest 
the presence of unshocked ejecta materials in the central region of the remnant (Rho et al. 
2008; Smith et al. 2009). However, it has been difficult to estimate how much amount of dust 
is there in the unshocked region from the observations covering the wavelengths of A < 70 
/xm, as pointed out by Rho et al. (2008). Since the mass of dust in Cas A is likely to be 
dominated by the unshocked low-temperature dust, the comprehensive far-IR observations at 
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A > 50 jjm are necessary to reveal the mass of dust formed in the Cas A SN. Q| Unfortunately, 
Cas A suffers from the contamination from the interstellar dust in the foreground molecular 
clouds, but the Herschel will be able to give more critical clues on the mass of cold dust in 
Cas A. 
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Fig. 1. — Relative number abundance of elements in the ejecta of the SN lib model at 600 
days after the explosion as a function of enclosed mass. The decay of radioactive elements is 
taken into account, and the original onion-like structure is assumed to be preserved without 
any mixing of elements. The progenitor mass of the SN lib is 18 Mq, and the mass at the 
explosion is 4.52 Mq. The mass cut is 1.58 M©, which results in the ejecta mass of 2.94 Mq. 
The mass of the hydrogen envelope is 0.08 Mq. 
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Fig. 2. — Structures of (a) temperature and (b) density of the gas within the He core at 
day 300 {solid lines) and day 600 {dashed lines) after the explosion, and (c) the velocity 
distribution for the SN lib model {thick lines). For comparison, shown are those for the SN 
II-P model {thin lines) with Meje = 17 and M^nv = 13.2 Mq (Umeda & Nomoto 2002). 
The mass coordinate is normalized by the mass of the He core. Note that both SN models 
have Ekin = 10^^ ergs and M{^^m) = 0.07 Mq. 
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Fig. 3. — (a) Mass distribution and (b) average radius of each dust species formed in the 
ejecta of the SN lib model as a function of enclosed mass. 
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time (yr) 

Fig. 4. — Time evolution of radius [upper panel) and velocity [lower panel) of the reverse 
shock for the uniform CSM models AO.l [dotted), Al [solid), and AlO [dashed) with nH,o = 
0.1, 1, and 10 cm~^, respectively. The trajectory of the boundary between the He core and 
the hydrogen envelope is depicted by the thin solid curve in the upper panel, where the 
boundary is assumed not to be decelerated after the collision with the reverse shock for 
illustration. 
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Fig. 5. — Structures of density {upper panel) and temperature {lower panel) of the gas at 
given times within the Type lib SNR for the model AlO with nu,o = 10 cm~'^. The positions 
of the forward and reverse shocks are indicated by the upward- and downward-pointing 
arrows in the lower panel, respectively. 
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Fig. 6. — Time evolutions of the positions [upper panel) and radii [lower panel) of dust 
grains within the Type II SNR for the model AlO. The initial radii of grains shown in the 
figures are 0.001 /im for AI2O3 and FeS [dashed lines) and 0.003 fim for C, MgSiOs, Mg2Si04, 
Si02, MgO, and Si [solid lines). The trajectories of the reverse and forward shocks as well 
as the position of the boundary between the He core and the hydrogen envelope are depicted 
by the thick curves in the upper panel. 
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Fig. 7. — Time evolutions of the total dust mass inside the Type lib SNRs running into 
(a) the uniform CSM with nH,o = 0.1 cm~^ (dotted), 1 cm~^ (solid), and 10 cm""^ (dashed), 
and (b) the stellar wind CSM with nH,i = 30 cm~^ (dotted), 120 cm~^ (solid), and 200 cm""^ 
(dashed). 
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Fig. 8. — Time evolution of the mass of each dust species inside the Type lib SNRs running 
into the stellar wind CSM with the density profile of p oc for (a) nn.i = 30 cm~^ and {b) 
nn 1 = 120 cm~'^. The thick solid line indicates the time evolution of the total dust mass. 
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Fig. 9. — Time evolutions of the gas density [upper panel) and the corresponding total ion 
number density [lower panel) as a function of distance from the center for (a) ?7,h,i = 30 cm""^ 
and [h) ?t,h,i = 120 cm~'^. The positions of the reverse shock are indicated by the arrows in 
the upper panels. In the lower panels, the positions of C [circles), Si02 [squares), and Si 
grains [triangles) are plotted where the filled and opened sympols denote the unshocked and 
shocked ones, respectively. 
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Fig. 10. — Time evolutions of the positions {upper panel) and radii {lower panel) of C, Si02, 
and Si grains with the initial radii of 0.003 /xm in the Type lib SNR for (a) nH,i = 30 cm""^ 
and (6) nH,i = 120 cm~^. The trajectories of the reverse and forward shocks as well as the 
position of the surface of the He core are depicted by the thick curves in the upper panel. 
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Fig. 11. — Time evolution of SEDs by thermal emission from shocked- heated dust within 
the Type lib SNR for the uniform CSM model Al (nn.o = 1 cm~^) every 200 yr after 
the explosion. The solid and dashed lines show the IR thermal emission with and without 
stochastic heating, respectively. The data points are the syncrotron-subtracted flux densities 
observed for Cas A and are taken from Hines et al. (2004). 
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Fig. 12. — Same as Figure 11, but for the model AlO (nH,o = 10 cm ^). 
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Fig. 13. — Same as Figure 11, but for the stellar wind CSM model B30 (?t.h,i = 30 cm ^). 
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Fig. 14. — Same as Figure 11, but for the model B120 (nH,i = 120 cm 
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Fig. 15. — Comparison between the IR observations of Cas A and the calculated SEDs by 
thermal emission from shocked dust at 330 yr for (a) nH,i = 30 cm~^, (6) 120 cm~^, and 
(c) 200 cm~^. The dashed lines indicate the contribution of each dust species. The IR flux 
densities in (a) are produced only by C grains. The data points are the same as in Figure 
11. 
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Fig. 16. — Contributions from (a) CS amorphous carbon and (6) CS astronomical silicate 
to the IR SED calculated for rin,! = 120 cm~^. The dashed and dotted lines denote the 
thermal emissions from the shocked CS dust and shocked newly formed dust, respectively, 
and the solid lines represent the total SED. The observation data for Cas A are the same as 
in Figure 11. 
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Fig. 17. — Contributions from the unshocked ejecta-dust to the IR SED calculated for 
TT-H,! = 120 cm~^ at (a) 330 yr and [b) 340 yr. The dashed and dotted lines denote the 
thermal emissions from the unshocked and shocked newly formed dust, respectively, and the 
solid lines represent the total SED. The temperature of cold dust is 43 K in (a) and 40 K in 
(b). The observation data for Cas A are the same as in Figure 11. 
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Table 1. Mass of the Main Refractory Species in the Ejecta of the SN lib 



species MA,i (Mq) M^^ (Mq) M^f /M^,, 



(1) 


(2) 




(3) 




(4) 


c 


1.14 X 10- 


-1 


7.08 X 10- 


-2 


0.621 





6.86 X 10- 


-1 


4.35 X 10" 


-2 


0.063 


Mg 


1.07 X 10- 


-1 


2.05 X 10- 


-2 


0.192 


Al 


9.31 X 10- 


-3 


3.28 X 10- 


-5 


0.004 


Si 


1.07 X 10- 


-1 


3.12 X 10- 


-2 


0.292 


S 


3.33 X 10- 


-2 


5.35 X 10- 


-4 


0.016 


Fe 


7.92 X 10- 


-2 


9.35 X 10" 


-4 


0.012 


others 


1.61 X 10- 


-1 










Total 


1.30 




0.167 







Note. — Col. (1): Species of the main heavy ele- 
ment in the ejecta of the SN lib. Col. (2): Mass of 
each heavy element (M^ j). Col. (3): Mass of each 
element that is locked up in dust grains at the time 
of dust formation (M^^^*). Col. (4): Mass fraction of 
each element that is locked up in dust grains at the 
time of dust formation {M'^^^ / M A,i) ■ 
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Table 2. Mass of Each Dust Species Formed in the Ejecta of the SN lib 



dust species 




) 






(1) 


(2) 




(3) 


(4) 


C 


7.08 X 10' 


^2 


0.423 


0.054 


AI2O3 


6.19 X 10' 


-5 


3.7 X 10"^ 


4.8 X 10-^ 


MgSiOs 


5.46 X 10" 


-2 


0.326 


0.042 


Mg2Si04 


1.74 X 10" 


-2 


0.104 


0.013 


Si02 


1.57 X 10- 


-2 


0.094 


0.012 


MgO 


2.36 X 10- 


-3 


0.014 


1.8 X 10-3 


FeS 


1.47 X 10" 


-3 


0.009 


1.1 X 10-3 


Si 


5.07 X 10" 


-3 


0.030 


3.9 X 10-3 


Total 


0.167 




1 


0.128 



Note. — Col. (1): Species of dust formed in the 
ejecta of the SN lib. Col. (2): Mass of each grain 
species (Mdj). Col. (3): Ratio of mass of each dust 
species to the total dust mass (M^j/M^). Col. (4): 
Ratio of mass of each dust species to the total metal 
mass (M^j/Ma). 
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Table 3. Results of the dust evolution calculation for each CSM model 



CSM model 
(1) 


j-He 
''coll 

(yr) 

(2) 


'-coll 

(yr) 

(3) 


''coll 

(yr) 

(4) 


-Mgui-v 

(Me) 
(5) 


V 

(M,,„/Md) 
(6) 


D330 
-"rev 

(pc) 
(7) 


0330 
Hovw 

(pc) 
(8) 


M330 

warm 

(Me) 
(9) 


(Me) 
(10) 


AO.l 


335 


1300 


3400 


1.3 X 10^^ 


7.8 X 10-^ 


2.76 


3.46 


0.0 


0.167 


Al 


155 


590 


1550 


0.0 


0.0 


1.98 


2.52 


1.0 X 10"^ 


0.167 


AlO 


75 


270 


700 


0.0 


0.0 


1.32 


1.83 


0.038 


0.094 


B30 


63 


340 


1850 


0.0 


0.0 


1.51 


2.20 


0.048 


0.111 


B120 


48 


150 


560 


0.0 


0.0 


0.93 


1.67 


0.008 


0.072 


B200 


45 


120 


370 


0.0 


0.0 


0.70 


1.50 


0.028 


0.011 



Note. — Col. (1): Name of the CSM models considered in this paper. The labels A and B 
represent the uniform CSM and the stellar wind CSM (p oc r~^), respectively, and the attached 
number denotes the value of nH,o for model A and nH,i model B (see text). Cols. (2), (3), and (4): 
Times at which the reverse shock encounters the He core, 0-rich layer, and Si-rich layer, respectively. 
Col. (5): Mass of surviving dust ai t = 10^ yr. Col. (6): Ratio of surviving dust mass Mgurv to the 
initial dust mass Md- Cols. (7) and (8): Radii of the reverse and forward shocks at t = 330 yr. Col. 
(9): Mass of warm dust heated up by the shock at t = 330 yr. Col. (10): Mass of unshocked cold 
dust at t = 330 yr. 



